SELF-AVERAGING SCALING LIMITS OF TWO-FREQUENCY WIGNER 
DISTRIBUTION FOR RANDOM PARAXIAL WAVES 



ALBERT C. FANNJIANG 

Abstract. Two- frequency Wigner distribution is introduced to capture the asymptotic behavior 
of the space-frequency correlation of paraxial waves in the radiative transfer hmits. The scahng 
hmits give rises to deterministic transport-hke equations. Depending on the ratio of the wavelength 
to the correlation length the limiting equation is either a Boltzmann-hke integral equation or a 
Fokker-Planck-like differential equation in the phase space. The solutions to these equations have 
a probabilistic representation which can be simulated by Monte Carlo method. When the medium 
fluctuates more rapidly in the longitudinal direction, the corresponding Fokker-Planck-like equation 
can be solved exactly. 



1. Introduction 

A central quantity of wave propagation in random media is the correlation of wave field at two 
space-time points. Through spectral decomposition of the time-dependent signal the space-time 
correlation is equivalent to the space-frequency correlation of the wave field [16] . The main focus 
of the present work is to derive rigorously closed form equations governing the space-frequency 
correlation of paraxial waves in the radiative transfer regime. 

For optical wave propagating through the turbulent atmosphere, the complex-valued wave am- 
plitude is governed by the stochastic Schrodinger (paraxial) equation with a white-noise potential 
[17j. The conventional approach uses the two- frequency mutual coherence function and various ad 
hoc approximations [13]. Recently, we introduced the two- frequency Wigner distribution in terms 
of which we derived rigorously a complete set of two-frequency all-order moment equations and 
solved exactly the mutual coherence function in the geometrical optics regime [7]. 

In this paper, we consider the different regime of radiative transfer and prove the self-averaging 
convergence of the two- frequency Wigner distribution for the paraxial wave equation. In other 
words, in radiative transfer, the whole hierarchy of two-frequency moment equations is reduced to 
a single radiative-transfer- like equation. 

Let Lz and Lx be, respectively, the characteristic macroscopic length scales of the wave beam 
in the longitudinal and transverse directions. We assume the phase speed in the vacuum is unity. 
The Fresnel number commonly defined as 



with the central wavelength Aq = 27: /k^ measures the significance of Fresnel diffraction. In all 
the scalings considered here the corresponding Fresnel number is large, indicating strong Fresnel 
diffraction effect. 

We use \o,Lz,Lx to non-dimensionalize the paraxial wave equation [13]. Let ki,k2 be two 
(relative) wavenumbers nondimensionalized by the central wavenumber fco- Then the wave field 
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of wavenumber kj satisfies 



where 7 is the reciprocal of Fresnel number, V represents the refractive index fluctuation with the 
correlation lengths iz and ix in the longitudinal and transverse directions, respectively, and fi is 
the magnitude of the fluctuation. We denote the ratios iz/Lz and ix/L^, respectively, by pz and 
Px which then are the correlation lengths of the medium fluctuations in the unit of characteristic 
length scales of the wave beam. 

Let us now describe a family of scaling limits which would yield non-trivial self-averaging limit for 
the space-frequency correlation of the wave beam. The main characteristic of the scaling limits is 
that the transverse correlation length ix is much smaller than the beam width Lx but is comparable 
or much larger than the central wavelength Aq. Roughly speaking, the wave beam experiences the 
spatial diversity as it propagates through the medium but it does not do so over a single wavelength 
in order to avoid the effective medium regime. This condition is best described by dimensionless 
quantities as 

(2) 9 ^ -f/px = 0(1) 

which includes the possibility that ^ ^ 1. This is the radiative transfer regime for monochromatic 
waves [8]. 

To reduce the complexity of technique, however, we restrict ourselves mainly to the longitudinal 
case where Pz/px = 0(1), including the possible scenario Pz/px *C 1. In this case, the medium 
fluctuation of the following order of magnitude 

(3) p = 



y/Pz y/Pz 

would lead to non-trivial multiple scattering effect (see [6], [8] and references therein). 

For wave fields of two different frequencies to interfere coherently, the unsealed frequency shift 
needs to be limited to 0(1) regardless the central frequency. In other words, we look for 0(1) 
coherence bandwidth which is a characteristic of the random medium. Again, this is best expressed 
in the dimensionless quantities as 

(4) lim ki = lim k2 = k, lim '~f^^k~^(k2 — ki) = 3. 

We assume /? > below. We call this the two-frequency radiative transfer scaling limits. 

Let us now review the basic framework of two-frequency Wigner distribution in the paraxial 
regime. The two-frequency Wigner distribution is defined as [7] 



where the scaling factor ^ykj is introduced so that W satisfies a closed- form equation (see below). 
Note here the scaling factor for the parabolic wave is different from that for the spherical wave 
introduced in |10j . 

The following property can be derived easily from the definition 

\\W{z,;-)h = (^^) ^ ||^l(^,-)l|2||*2(^,-)||2. 
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Since "ifj are governed by the Schrodinger-like equation the L^-norm does not change with z, i.e., 
||VF(2:, •, •)||2 = ||VK(0, •, •)||2. The Wigner distribution has the following obvious properties: 

7y 
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n'^Vkik2\'^ , p\/k\ 



+ 



) 



7 y -47 2 ' 47 2 

and so contains essentially all the information in the two-point two-frequency function. 
The Wigner distribution Wz satisfies the Wigner-Moyal equation exactly 



(8) 



dW 1 

+ p • Vxl^ + ^^^W^ = 



d 



where the operator Lz is formally given as 



CM 



3-1 



e'^^-^l^^kxW^z, X, p + -^) - e^'i-*/^A;2VF(z, x, p - 

-1/2, 



Pz 



with X = x/^Oj; being the 'fast' transverse variable. As pz 0, pz Cz displays the classical central 
limit scaling in the z-variable and thus this is referred to as the longitudinal case. The complex 
conjugate VF^*(x, p) satisfies a similar equation 



(9) 

where 



dW* 1 

+ p • VxVK* + ^C*zW* = 



dz 



3-1 



ei^-^/Vk^k2W*{z,^,p + 



e'n-s^/VkIkiW*{z,^,p 



2Vk2 

We use the following definition of the Fourier transform and inversion: 
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y(-,dq). 

Pz 



{2ttY 



-jxp 



/(x)dx 



g{p)dp. 



When making a partial (inverse) Fourier transform on a phase-space function we will write J-'i (resp. 
J-i^) and J-2 (resp. ^2^) ^° denote the (resp. inverse) transform w.r.t. x and p respectively. 

2. Formulation and theorems 

To describe the scaling limits more efficiently we introduce two controlling parameters e,a > 
and set 

(10) pz = e^, Px = e'". 

By ©-(U]) all the other parameters can be expressed in terms of e, a and 9. The radiative transfer 
scalings in the longitudinal case then correspond to 



(11) 



e ^ 0, a G (0,11, lim 6* < 00 

e— >0 



along with ©-(g]). 

To get the self- averaging result, it is essential to consider the weak formulation of the Wigner- 
Moyal equation: To find Wz in the space C([0, cxd); L^(R^'^) of z-continuous, L^(IR^'^)-valued pro- 
cesses such that llM/'f II2 < llM^lb^Vz > 0, and 

(12) {w^,@)-{Wo,e) = r {w^,p-v^Q)ds + - r {w^,c:e)ds, 

Jo £ Jo 



where Wq £ L'^{M?'^) is the initial condition and we consider all the smooth test functions of compact 
supports G C^(M^'^). Here and below L^(M^'^) is the space of complex- valued square integrable 
functions on the phase space M^*^ endowed with the weak topology and the inner product 

{Wi,W2) = J H^i*(x,p)M^2(x,p)dxdp. 

We define for every realization of the operator £* to act on a phase-space test function Q as 

(13) /::e(x,p) ^ -i9-'j^2 [<5.v;^(x,y)^2~'0(x,y)] 

with the difference operator 5^ given by 

(14) SeV:{^,y) = feiy/(^ + -^)-fc2y/(-'' ^ 



We define Cz in the similar way. 

The existence of solutions in the space C([0, oo); L^(M^'^)) can be established by the same weak 
compactness argument as in [7j. We will not, however, address the uniqueness of solution for the 
Wigner-Moyal equation ()12p but we will show that as e — > any sequence of weak solutions to eq. 
()12p converges in a suitable sense to the unique solution of a deterministic equation (see Theorem 1 
and 2). 

We assume that 14 (x) = V{z,x) is a real- valued, centered, z-stationary, x-homogeneous ergodic 
random field admitting the spectral representation 



Vz{x) = j exp(fp-x)V;((ip) 



with the z-stationary spectral measure V^(-) satisfying 

E[T4(dp)T4(dq)] = <^(p + q)$o(p)dpdq. 

The transverse power spectrum density is related to the full power spectrum density $(t/;,p) as 
*^'o(p) = / ^{w, p)dw. The power spectral density $(k) satisfies ^>(k) = $(-k), Vk = {w, p) € M'^^^ 
because the electric susceptibility field is assumed to be real-valued. Hence ^{w,p) = <&(— w,p) = 
^{w, — p) = ^{—w, — p) which is related to the detailed balance of the limiting scattering operators 
described below. 

The first main assumption is the Gaussian property of the random potential. 

Assumption 1. y(z,x) is a Gaussian process with a spectral density <I>(k),k = {w,p) G M'^"^^ 
which is smooth, uniformly bounded and decays at \k.\ = oo with sufficiently high power of \k.\~^ . 

We note that the assumption of Gaussianity is not essential and is made here to simplify the 
presentation. Its main use is in deriving the estimates ()40p . (I4ip below. 

Let J^z and be the sigma-algebras generated by {Fg : Vs < z} and {Vs '-^s > z}, respectively 
and let L^{J^z) and L^{J^^) denote the square-integrable functions measurable w.r.t. to them 
respectively. The maximal correlation coefficient p{t) is given by 

(15) p{t) = sup sup E [kg] . 

E[h]=0,E[h2] = l j.[^,^o,Eb2] = l 

For Gaussian processes the correlation coefficient p{t) equals the linear correlation coefficient 
given by 



(16) sup / R{t - Ti - T2,'k)gi{Ti,'k)g2{T2,'k)d'kdTidT2 

91,92 J 



where R{t,k) = J e^*?$(C,k)d^ and the supremum is taken over aU gi,g2 G which are 

supported on (—00, 0] x and satisfy the constraint 

R{t-t','k)gi{t,k)gl{t',k)dtdt'dk= j R{t - t' ,k)g2{t,k)g^{t' ,k)dtdt' dk = 1. 



There are various criteria for the decay rate of the hnear correlation coefficients, see |12] . 
Next we make assumption on the mixing property of the random potentiaL 

Assumption 2. The maximal correlation coefficient pit) is integrable: p{s)ds < 00. 

We have two main theorems depending on whether hm^^o is positive or not. 

Theorem 1. Let 9 > be fixed. Then under the two-frequency radiative transfer scaling 
U0\} the weak solutions, denoted by W^, of the Wigner-Moyal equation (0j converges in probabil- 
ity in C([0, 00), L^(]R'^)), the space of L"^ -valued z-continuous processes, to that of the following 
deterministic equation 

(17) ^W,+p.VW, = ^lK{p,q) e~*^^'i-''/(2v^)M^,(x,p + -^)-t^,(x,p) dq 

where the kernel K is given by 

K(p,q) = «3>(0,q), for qG(0,1), 

and 

i^(p,q) = $((p + ^) •q,q), for a = 1. 

Theorem 2. Assume hm^^^o ^ = 0. Then under the two-frequency radiative transfer scaling f^- 
0; ( fiOj) the weak solutions of the Wigner-Moyal equation (0j converges in probability in the space 
C([0, 00), -L^(M'^)) to that of the following deterministic equation 

(18) ^W, + p-VW, = k(Vp- ^/3x^ • D • (^Vp - ^/3x^ 
where the (momentum) diffusion coefficient D is given by 

(19) D = 7ry^>(0,q)q(g)qdq, for a & (0,1), 

(20) D(p) = Try $(p • q, q)q(g) qdq, for a = 1. 

Remark 1. In the transverse case of px <^ Pz (or a > 1), then with the choice of p = ^fp^ (or e"^) 
the limiting kernel and diffusion coefficient become, respectively 

(21) i^(p,q) = S{ip + ^)-q) I ^w,q)dw, 

D(p) = vtIpI^"'" / / ^{w,p±)dw p± p±dp±. 

J p px=0 J 

The proof of such result requires additional assumptions which would complicate our presentation, 
so we will not pursue it here, cf. [8], [6]. 



Remark 2. The form of ( f j7| ) and 1^1^ suggests the new quantity 

2B,(x,p) =e-¥-PTF,(x,p) 



in terms of which the equations can be written respectively as 



2IT,(x,p + ;^)-21J,(x,p) 



dq, 



(22) A2B, + p.V2H, + ||pp2n. = ^|i^(p,q) 

(23) |-2IJ, + p.V2ir, + ^|p|2!m, = fcVp-D-Vp2B,. 

oz 2 

The advantage of this is that the solution then is amenable to probabilistic representation as the 
scattering terms on the right side of equations are non-positive definite. Let {x{z),p{z)) be the 
stochastic process with x(0) = x, p(0) = p generated by the operator — p • Vx + A where A is the 
operator on the right hand side of either i22\) or 123\) . Then by Dynkin's formula we have 

(24) 2B,(x,p) =Ex,p [e-^/o IPWI''^*e-¥''(^)-P(^)H^o(x(z),p(z)) 

where Wq is the initial data for the two-frequency Wigner distribution and Ex,p is the expectation 
with respect to the probability measure associated with {x{z),p{z)). The probabilistic solution \24^ 
can be simulated by Monte Carlo method. 

When /ci = A;2 or /? = 0, eq. (fT7|) and (fTSj) reduce to the standard radiative transfer equations 
derived in [6], [8]. In view of the definition ([5]) the two-frequency Wigner distribution, captures the 
space- frequency correlation on the nondimensionahzed scale of 7 (in space as well as in frequency). 
The self- averaging property simply reflects the fact that on the wave field essentially decorrelates 
on the macroscopic scale ^ 7 and the rapidly fluctuating limit is then statistically stablized by 
coupling with a smooth deterministic test function. 

A notable fact is that eq. (jlSp with (jl9p is the same governing equation, except for a constant 
damping term, for the ensemble- averaged two-frequency Wigner distribution for the 2;-white-noise 
potential in the geometrical optics regime, [7]. Fortunately eq. (jl8p - (|19p is exactly solvable and the 
solution yields asymptotically precise information about the cross-frequency correlation, important 
for analyzing the information transfer and time reversal with broadband signals in the channel 
described by the stochastic Schrodinger equation with a z- white- noise potential [9] (see also [2], [3], 
[3]). Moreover, eq. (|18p -(|19 p arises as the boundary layer equation for the two-frequency radiative 
transfer of spherical wave [lOj. 

The proof of Theorem 2 is analogous to that of Theorem 1 and for the sake of space we will 
not repeat the argument. We refer the reader to [6] for the needed minor modification. More 
directly eq. (jlSp can be obtained from eq. (jl7p in the limit — > by Taylor expanding the terms 
involving 9 on the right hand side of (jl7p up to second order in 9. The first-order- in-0 terms are 
also first-order-in-q and thus vanish due to the symmetry i^(p, q) = i^(p,— q). The remaining 
terms, after dividing by 9^ and passing to the limit, become the right hand side of (llSp . 

For the proof of Theorem 1 below, we set 9 = \ for ease of notation. 



3. Martingale formulation 

For tightness as well as identification of the limit, the following infinitesimal operator will 
play an important role. Let = V{z/e'^,-). Let J^l be the cr-algebras generated by {V7, s < z} 
and E| the corresponding conditional expectation w.r.t. T^. Let be the space of measurable 
function adapted to {.F|,Vt} such that sup^^-^g E|/(z)| < co. We say /(•) S V{A^), the domain of 
A^, and A^f = gii f,gGM^ and for f\z) = (5-i[E|/(z + 5)- f{z)] we have 

supE|/(z)| < 00 

z,5 

limE\f^{z)-g{z)\ = 0, Vz. 
0— >o 
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Consider a special class of admissible functions f{z) = cj}{{W^ , Q)) , f (z) = </)'((VFJ, B)), V(/) G 
C°°(M) we have the following expression from (jl2p and the chain rule 



(25) 



f'iz) 



{W!,p-VQ) + -^{W!,£:@) 



In case of the test function that is also a functional of the media we have 



(26) 



A'fiz) 



(H^f ,p • ve) + i {w!,i:iQ) + {w!,A'e) 



and when Q depends explicitly on the fast spatial variable 
the gradient V is a sum of two terms: 



V = Vx + e" 



-2a 



Vx 



where Vx is the gradient w.r.t. the slow variable x and Vx the gradient w.r.t. the fast variable x. 
A main property of A"^ is that 



(27) 

Also, 

(28) 



f{z) - / A'f{s)ds is a J^|-martingale, V/ G V{A' 



WJ{z)- f{s) = j\lA'f{T)dT ^s<z a.s. 



(see [H]). We denote by A the infinitesimal operator corresponding to the unsealed process Vz{-) 
V{z,-). 



4. Proof of tightness 
In the sequel we will adopt the following notation 

/(z)^0((H^,^e)), nz)^cp\{w!,e)), r(z)^0"((i^|,e)), 



Namely, the prime stands for the differentiation w.r.t. the original argument (not t) of /, /' etc. 

A family of processes {W|,0 < e < 1} in the Skorohod space /^([O, oo); L^(M^'^)) is tight if and 
only if the family of processes {{W^,e),0 < e < 1} C D{[0,oo); Ll(M?'^)) is tight for all 6 G 
[llj . We use the tightness criterion of [15j (Chap. 3, Theorem 4), namely, we will prove: Firstly, 



lim limsupP{sup | {W^, G) | > iV} = 0, Vzq < oo. 

N^oo z<zo 



(29) 



Secondly, for each <f> G C°^(M) there is a sequence f^{z) G V^A^) such that for each zq < oo 
{A"^ /'^{z), 0<e<l,0<2;< zq} is uniformly integrable and 



(30) 



limP{sup |/=(z) -0((Wf,e))| >5} = 0, V5>0. 

z<zo 



Then it follows that the laws of G) , < e < 1} are tight in the space of -D([0, oo); M). 

To prove the tightness in the space C([0, oo); ^^(M^'^)) let us recaU that G C([0, oo); LliM?"^)) 
and that the Skorohod metric and the uniform metric induce the same topology on C([0, oo); L'^{M?'^) 

Let us note first that condition (j29p is satisfied because the L^-norm is uniformly bounded. 
The rest of the argument for tightness will be concerned with establishing the second part of the 
criterion. 
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Consider now the expression 
(31)£:e(x,x,p) = -ie~^ 



A;ie(x, p 



e^'i-^/^fc2e(x,p + 



q 



or equivalently 

where SeVg is defined by ()14p . It is straightforward to check that (j3ip solves the corrector equation 
(33) [e^^ap . y _ ^ ^ej ^ e-2£*e 

Recall that Vx and Vx are the gradients w.r.t. the fast variable x and the slow variable x, 
respectively. 

We have the fowllowing estimate. 

Lemma 1. 



E 



(x,p) 



< 



r p{s)ds] [ e*'i-*/^A:ie(x, p + - e^'i-^/^A:2e(x, p 

Jo J J 2yjki 



'$(e,q)(iC(iq. 



Proof. Consider the following trial functions in the definition of the maximal correlation coefficient 
h = /is(x,p) 

„iq.x/v/fcr^^0(x,p + _ e^^-*/^fc2e(x,p e''="'(^-")P-'i/"'"E^l//(dq) 



2Jh 



2Vko 



9 = 9t{x,P) 

= i 



2\/ki 2\/k2 -I 

It is easy to see that hg G L'^{P,^},J^^-2^)gt GG L'^{P,Q,,J-^_2^) and their second moments are 
uniformly bounded in x, p, e since 

E[/i2](x,p) < E[<72](x,p) 



E 



b,2](x,p) = / e^'i-/^A;iG(x,p + -^)-e^'i-/v^A;2G(x,p--^)'$(e,q)dCdq 



which is uniformly bounded for any integrable spectral density <I>. 
From the definition (1151) we have 



|E[/i,(x,p)/i,(y,q)]| = |E [hM < p{e-\t - z))W}l^ [/i^(x,p)] W}/^ [5t'(y,q)] ■ 
Hence by setting s = t,x = y,p = q first and the Cauchy-Schwartz inequality we have 

E[/i2(x,p)] < p\e-\s - z))ngl{-^,v)] 
|E[/i,(x,p)/it(y,q)]| < p{e-\t- z))p{e-\s- z))¥}%l{-K,p)\¥}%l{y,^)\, Vs,t>^,Vx,y. 
Hence 



E[/i,(x,p)5t(x,p)]dsdt < E[52](x,p) 
which together with p4p yields the lemma. 



p{s)ds 



□ 



Corollary 1. 











" poo 




(34) 


E 


P • Vx£:e 


'(x,p) < 


1 p{s)ds 
Jo 
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e*q-x/v^^^p.v^e(x,p + -^; 

-e'^-^/^^k2P ■ Vxe(x, p - q)d^dq. 



Inequality (I34p can be obtained from the expression 
p • Vx£*0(x,x,p) = ie~ 



*^-/v^A;2p-Vxe(x,p + -^ 



2a/ A; 



as in Lemma [TJ 

We will need to estimate the iteration of and CI: 



£:£:G(x,i,p) 



1 /"OO 



giq-x/VfeTgiq'-x/Vfci 



x/cfe(x,p + + 



q 



q ^ q 



2^/k2 2^x' 



,j(s-2)q'-q/(2e2'^) 



xe 



■g.qx/v/teg.q'.x/V^^^^^0(^^p^ q _ 
-i(s-2)q'.q/{2e2")\ 



q ) _ e*i*/v^e*'i'-'^/^A;iG(x,p - 



q 



/:*/:*e(x,x,p) 

/>oo />oo 



dsdt I E^V7((iq)E^[l//(dq')]e^(^"")P-'i/" \iit~^)p-<i/^ 

kik2<d{x, p 



gjq-x/v^giq'-x/v^ 



+ 



2Vfc^ 2y/h' 



,j(s-2)q'-q/(2e2<^) 



q 



xe 



giqx/Vfc^giq'.x/v^^^^^Q^^^ p ^ _q q ) _ g*qx/v^giq'.x/v^^20^^^ p _ _q 

2^/kl 2\/k2 2\jk2 2\jk2' 

-i(s-^)q'.q/{2e2-)1 



which can be more easily estimated by using ([32]) as follows. First we have the expressions after 
the inverse Fourier transform 

(35)^-2-1 {£:£:e}(x,i,y) = e'^ y 5,l//e— ""(^-^)^--^'' [E,[<5,y/].F2-iG] (x, x, y)ds 



(36)^^2-1 {£:£:e}(x,i,y) 



_ie-2a(t_^)Vy.Vs 



[E,[d,y/].F2-ie]}(x,i,y)(is(it. 



Lemma 2. 



-e*'i-/v^A;2e(x,p--^; 



^{£,,q)dS,dxdqdp 



for some constant C independent of e. 

Proof. Let us consider £*£*0 here. The calculation for C*C*Q is similar. 
By the Parseval theorem and the unitarity of exp (irVy • Vx), r E M, 



< 



Cos-' J |E{E,[5,F/]e— "'"(^-^)^--^^ [E,[5,Vn:F,'e]] i^,5c,y) 



dsdt 



E 



ds'dt'dxdy 



+Coe-'l |E{e— -'"(^-^)^--^- [E,[d,y/].F2"ie] (x,y) 



xe 



^'-)v,.v, [^^ls^ye^^Jr^lQ-^ (x,i,y)} 



dsdtds'dt'dxdy. 



The last inequality follows from the Gaussian property. Note that in the x integrals above the fast 
variable x is integrated and is not treated as independent of x. 
Let 



9it) = SeVi 



and 



h{s) = e-'' '"(^-^)VyVs [5eV!r^^Q] . 
The same argument as that for Lemma [1] yields 

|E[E,[5(t)]E,[/i(s)]]| < EV2[E,[5(t)]2]Ei/2[E,[/j(,)]2] 

< p{E-\t-z))p{e~\s-z))¥}'''[g\t)]¥}/^[h\s)l t,s>z; 
\E[EMmz[9{t')]]\ < EV2[E,[5(t)]2]Ei/2[E,[5(t')]2] 

< p{e~\t-z))p{e-\t'-z))E^/'[g\t)]E^/^[g\t% t,t' > z; 

\E[EMs)mhis')]]\ < E'/^[EMs)?]^'^^mHs')?] 

< p{£~\s- z))p{£~\s' - z))E^/^[h'^{s)]E^/^[h\s')], s,s'>z. 
Combining the above estimates we get 

r r „ . . -,2 

dxdy 



E\\cic*Ml < p{s)dsj y e[4f/]2e [e— ~ "(^-^)^--^^ [5,y/.F2-ie] 

< C2 (^1^^ p{s)dsy E[Vn' I |e^'i-''/^A:iVxe(x,p + ^) 



-e^'i-*/^fc2Vxe(x, p ^) $(^, q)dCdxdcidp 

2v^2 
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□ 



Now let us consider the second moment of p • VyiC*,C*,Q and 
^2-i{p-Vx£:£:G}(x,i,y) 



26 



The same calculation as in Lemma [2] yields the following estimates: 
Corollary 2. 

2\/fe2 



-E[y/]2 / e^'i-Zv^fcip • Vxe(x, p + 



E||£:£:£:eii^ - ^' /o ^'^^'^ 



ei4 



2VA;2 

jq-x/Vfcl 
2 



e"*"' ^■"^fci0(x, p + 



q 



2^Jk^ 



-e^^-*/^A;2G(x, p - $(^, qjd^dxdqdp 

2vfe2 



/or some constant C independent of e. 
Let 

(37) h{z)=ef'{z)(wlCie 
be the 1-st perturbation of f{z). 
Proposition 1. 

lim sup E|/i(z)| = 0, lim sup = in probability 



z<zo 



z<zo 



Proof. We have 

(38) 

and 

(39) 



E[|/i(.)|]<e||/'||oo||TVo||2E||£:e||2 



SUp|/f(z)| <e||/'||oo||T^0||2SUp||£:e||2. 
z<zo z<zo 
11 



Since C^Q is a Gaussian process and £*£*G is a x^-process, by an application of Borell's inequality 
PP we have 

(40) E[sup ||£:e||2] < C71og(-)E||£:e||i; 

z<zo e 

(41) E[sup ii£:£:ei|2] < ciogH-)E\\c:c:e\\i 

z<zo £ 

i.e. the supremum over z < zq inside the expectation can be over-estimated by a log(l/e) factor 
for excursion on the scale of any power of Hence the right side of (|38p is 0(e) while the right 
side of ([39j) is o(l) in probability by Chebyshev's inequality. □ 

Set f'^{z) = f{z) + fi{z). Then (j30p follows immediately from Proposition 1. 

Let us now prove the uniform integrability of A^f^. A straightforward calculation yields 

A'h = efiz) Iwi, p • Vx£:e\ + er\z) {W^ p • VxB) (w^ C*Q 



+nz) (wic*Aq) + f"{z) {wi,ciQ) (w^ciq) - -J'{z) {wixiQ) 

and, hence 

(42) A'fiz) = /'(z)(t^,^p.Vxe) + /'(z)(^^,^£:£:G) + r(^)(l^,^£:G)(vFf,£:e 



= Ao{z) + Ai{z) + A2{z) + Ri{z) 
where Ai{z) and A2{z) are the 0(1) statistical coupling terms. 
Proposition 2. 



lim supE\Ri{z)f = 

£~*Oz<zo 



Proof. First we note that 

\Ri\ < e 
Clearly we have 



^"lloollw^ollillP • Vxe||2||£:e||2 + Wfiumup ■ Vx(£:e)|h 



lim supE|i?i(z)|2 = 0. 

i^-^O z<zo 

by Lemma [J and Corollary [TJ □ 

For the tightness criterion stated in the beginnings of the section, it remains to show 
Proposition 3. {^^/^} are uniformly integrable. 

Proof. Let us show first that {^i}, i = 0, 1, 2, 3 are uniformly integrable. 
For this we have the following estimates: 

|^o(^)l < ll/'lloollvTolbllP-Vxeib 
\Aiiz)\ < \\f'\\^\\Woh\\c:cie\\2 
\A2{z)\ < ||/"||oo||H^ollil|/::e||2||£:G||2. 

The second moments of the right hand side of the above expressions are uniformly bounded as e ^ 
by Lemmas [1] and [2] and hence Aq{z), Ai{z), A2{z) are uniformly integrable. By Proposition [21 Ri 
is uniformly integrable. Therefore A'^f^ is uniformly integrable by (j42p . □ 
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5. Identification of the limit 



The tightness just estabhshed permits passing to the weak hmit. Our strategy for identifying the 
Umit is to show directly that in passing to the weak hmit the hmiting process solves the martingale 
problem with null quadratic variation. This would imply the limiting equation is deterministic. 
The uniqueness of solution to the limiting deterministic equation for given data then identifies the 
limit. 

For this purpose, we introduce the next perturbations /2, /a- Let 



(43) 

(44) 
where 

and 

Clearly, 

Let 
and 
Let 



4^^(V) = y"v'(x,p)Qi(G0e)(x,p,y,q)V(y,q)dxdpdyciq 
4'\V') = y"Q'ie(x,p)V'(x,p) dxdp, VV-GL^ 



Qi(e®G)(x,p,y,q) 



E 



/:*e(x,p)Z:G(y,q) 



Q;e(x,p) = E[£:£*e(x,p) 

Q2{@ ® e)(x, p, y, q) = E [£*M^, p)jC,l@{y, q) 
Q'2e(x,p)=E[£:£:e(x,p)" 

A^^\'iP) = y"v'(x,p)Q2(e®G)(x,p,y,q)V'(y,q) dxcZpcZydq 



Define 

(45) 

(46) 

Proposition 4. 



f2{z) 



-f"iz) 



Hz) = '-nz)[(wi,CAQ)-Af\wi) 



lim sup E|/2(2;)| = 0, lim sup 11/3(2;) | = 0. 



£^0 



Z<Zo 



z<zo 



Proof. We have the bounds 



aA/|) 



supE|/2(z)| < supe^liriloo \\Wo\\in^*M\i + nA'>{W- 

z<zo z<zo 

supE|/|(z)| < supe^ll/ll^ [||Wo||2E||£*£*e||2 + E[4')(W|)] 



z<zo 



z<zo 



A straightforward calculation shows that ¥,[A2 (W^)] and E[A'[^''(W|)] stay uniformly bounded 
w.r.t. e. The right side of the above expressions then tends to zero as e ^ by Lemma 1 and 
2. □ 
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l(2)/ 



We have 



A'hiz) = f"{z)[-{w^,c:e){w!,cie) + A^^\wi)]+R2iz) 

A'fsiz) = fiz)\-(w!,C:iC:Q))+A^^\w!)]+Rsiz) 



with 



R2(Z) 



2 

-2 f" I 



(W^|,p-VxG) + -(VFf,£:e) 



(47) 



w^,^p•Vx(/::e)) + -(H^,^£:/::e 



(2) 

where Gq denotes the operator 

CgV^ y" Q2(e0G)(x,p,y,q)V'(y,q)fiydq. 



Similarly 

(48) i?3(^) 



-eV'(^) 



w^l, p • Vx(/::/::e) ) + - ( wi, cicic^q 



(VFf,p-VxG) + -(VFf,£:G) 



VFf,/::£:G)-4')(Ty; 



(W^I,P • Vx(Q'2G)) + i (W^f ,£:Q'2G) 



Proposition 5. 



lim sup ^R2{z)\ = 0, lim sup E|i?3(z)| = 



2<Z0 



£^0 



Proof. Part of the argument is analogous to that given for Propositions) The additional estimates 
that we need to consider are the following. 
In R2: First we have 



sup £ E 

z<zo 



£2 / E[t^|(x,p)t^f(y,q)]E p • Vx£:G(x, p)£:G(y, q 



dxdydpdq 



E \Wl{^, p)M^f (y, q)] eV2[p . Vx£:G]2(x, p)EV2[£*e]2(y, q)dxdydpdq 



which is O(e^) by using Lemma [H Corollary [T] and the fact E [pyj(x, p)VFf (y, q)] G L^(M'^'^) in 
conjunction with the same argument as in proof of Lemma 1; Secondly, we have 



sup eE 

Z<ZQ 



supe||W^o||2E||£:E 

Z<ZO 



supe||VFo||2E||.F2-^£:E 

2<Z0 



£:G Ll<d 



Define 
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We then have 



= E 
< Ei/2 



j j J,F/(x,y)E [E,[/i,(x,y)]E,[/it(dx',dy')]] J^^^W^i^ ,y')d^ dy' dsdt 
%Vl{^,y)\p{e-\s - z))p{e-\t - z))¥}/^[hs{^,y)f 



j E[^t(dx',dy')]^dx'(iy'^ (^j (x', pOpdx'dp'^ dsdt 



dxdy > . 



Recah that ||W|||2 < ||Wo||2 and 

j E[/i4((ix', dy')fd^dy' = j [e(x, p + q/2) - e(x, p - q/2)]2$(^, q)d?dqdx(ip < oo 



so that 



e||^2-1£:e 

< llT^olbE^/'ll/^.i ( supE[5,y/]2 j y ^ p{e-\s - z))p{e-\t - z)) 

xp{e-^{s' - z))p{e-^{t' - z))E^/^\\hs\\2^'^^^\\hs'\\ldsdtds'dt' 



< ||Wo||2E3/2||/j,||2(supE[5,y/]2 



p{s)ds 



< oo. 



Recall from ([Ml) that 



E\\h "2 



Hence 



J [e(x, p + q/2) - e(x, p - q/2)]2$(e, q)dedqdxdp < 

= 0(e). 



CXD. 



sup eE 

z<zo 



In 



sup eE 

Z<20 



< e||T^o||2SupE||/::44G||2 



z<zo 



which is 0(e) by Corollary [2j 

The other two terms in have the respective expressions 

e2E|(H^,^p.Vx(Q'2e)>| < e^woh^'/^p-v^nc:c:Q]\\i 

which is O(e^) by Corollary [2] and 

eE\{W!,£:Q'2e)\ < e||W^o||2lE||/::E[£:£:e]||2 



< e||T^o||2(supEi/2|5,y/|2)Ei/2||4£*e|| 



which is 0{e) by Lemma [2j 



□ 



15 



Consider the test function f^{z) = f{z) + fi{z) + f2{z) + fz{z). We have 



(49) A^nz) 

= f'{z) {W!,p- V.e) + f"{z)A^^\w!) + f'A?{Wl) + R,{z) + R2{z) + Rs{z). 

Set 



(50) 



R^{z) = Ri{z) + R2{z) + Rsiz) 



It follows from Propositions 3 and 5 that 



Proposition 6. 



Proof. We have 



lim supE\R%z)\ = 0. 

£-»0 2<zo 



lim sup sup ylg^^ (-0) = 0, 



£^0 



Z<ZO \\M\2 = 1 



J V'(x, p)Qi(e (g) e)(x, p, y, q)^(y, q) dxdp dydq 
^ J V'(x, p)Qi(x,p,y,q)V'(y,q)dxdpdydq 



where Qi is defined by 



Qi (x, p, y , q) = [Qi (6 ® 6) (y , q, X, p) + Qi (e ® e) (x, p, y , q)] . 



The symmetrized kernel has the following expressions 



Qi(x,p,y,q) 

ds I dp'6(s,p')e*P'-(^-y)/^"'e-'''P-P'^' 



5'p'-*/^/ciVxe(x,p + -^) 



.gip'-x/Vfe^;.2Vxe(x,p--^; 



e'P'-*/^fciVxe(x, q + -^) - e^P'-/^^fc2Vxe(x, q 

2\/ki 



= 27r 



e^P'-"/^A;iVxe(x, p + -^) - e^P'-/^fc2VxG(x, p - -^) 

2^/kl 2\JK2 



e'P'-*/^A;iVxe(x, q + -^) - e^P'-*/^A;2Vxe(x, q - -^)1 $(p • p'e^-a^^ p')rfp'. 

2vA;i 2vA;2 -I 



which, as the inverse Fourier transform tends to zero uniformly outside any neighborhood of x = 
because of Assumption 1, and stays uniformly bounded everywhere. Therefore the L^-norm of 
tends to zero and the proposition follows. 
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Similar calculation leads to the following expression: For any real-valued, L^-weakly convergent 
sequence ip, we have 



lim AY'ii;' 



lim / ds I dwdcidxdp V^(x, p)$(?i;, q)e*''"e-*"P-'i^ 



ia-xe'''^"(— —) / 11 1 \ 

k,k2e^ 'VH v^^G(x,p + -(^ + ^)q) -A;2e(x,p) 



_|_g 



-iq-xe-2"(^ ^) / 1,1 1 

A;iA;2e v''-2 9 x, p 



fc^lim / ds / dwdqdxdp ^'^(x, p)$(it;, q)e*""'e-**P-'i^ 



^ 3|q| ^2-2q 

+e 2^ 



,-iq-x, 



/^/(2'=^^^)e(x,p--|)-G(x,p)]} 



where we have used (jH). Note that the integrand is invariant under the change of variables: s 
—s, q — > — q. Thus we can write 

lim A^^\i;') 

£— >0 



/tMim / ds I dwdqdxdp ij%x,p)^{w,q)e'''"e-''P-'^'' '° 



■^.q.x/3/(2fcV2)Q/ q 



e(x,p) 



27rA:2lim / dqdxdp V'(x, p)$(e2-2«(p + . q, q) re^^-^^/(2/cV2)Q^^^ p ^ q ^J _9(x,p) 



from which we obtain 
Alii;) = limA5^^(V'^) 



2Tre J dqdxdp V(x, p)$(0, q) [e''i-''^/(2'='^')e(x, p + ^) - e(x, p)j , a G (0, 1) 



2TTe f dqdxdp v(x, p)$((p + ^) • q, q) e^'i-^^/(2'='^')e(x, P + - e(x, p 



a = 1 



Recall that 



M|(e) = /(z) + /i(z) + /2(z) + /3(z)- r/'(^)(VF|,p-Vxe)ds 

JO 



(51) 



f"{s)A^^\w!) + f'{s)A^^\w!) ds- I R^s)ds 



is a martingale. The martingale property implies that for any finite sequence < zi < Z2 < < 
... < Zn < z, C^-function / and bounded continuous function h with compact support, we have 

(52) E{h{{w!^,e),{w!^,Q),...,{w!^,e))[Mi^,{e)-Mi{e)]} = o, 

Vs > 0, Zi < Z2 < ■ ■ ■ < Zn < Z. 



Let 



Af{z) = f'iz) [{W,,p-V^e)+Ai{W,)] . 
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Here and below, by slight abuse of notation, f{z) and f'{z) stand for (j){{Wz,Q)) and (p' {{Wz,Q)), 
respectively. In view of the results of Propositions 11121 [3l [H [H [6] we see that 

E{h{{w!^,Q),{w!^,e),...,{w!^,e))[r{z)-cp{{w!,e))]} = o, 
E{h{{w!^,e),{w!^,e),...,{wi,e))[A'f{z)-A^{{w!,e))]} = o. 

With this and the tightness of WJ we can pass to the limit e ^ in ([52|) . cf. [5j, Chapter 4, 
Theorem 8.10. Consequently that the limiting process satisfies the martingale property that 

E {h {{Wz„e) , {Wz„e) , e)) [M,+,(e) - M,(e)]} = o, > o. 

where 

M,(e) = /(z)- r Af{s)ds. 
Jo 

Then it follows that 

E[Mz+s{Q)-Mz{e)\Wu,u<z]=0, Vz,s>0 

which proves that Mz{Q) is a martingale. 
Choosing (j){r) = r and we see that 

M«(e) = (M^„e) - r [(Ty„p-Vxe) + ii(TF,)] ds 

Jo 

is a martingale with the null quadratic variation 



MW(e),MW(e) 



0. 



Thus 



Since (VFJ, 0) is uniformly bounded 

\{w!,e)\<\\Woh\m, 

we have the convergence of the second moment 

YimE[{w!,ef] = {Wz,ef 



e 

and hence the convergence in probability. 
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